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ABSTRACT 

Infrared Dark Clouds (IRDCs) host the initial conditions under which massive stars 
and stellar clusters form. It is therefore important to study the kinematics, as well as 
the physical and chemical properties of these regions. Their complex structure however 
posits challenges in the data interpretation. We have obtained high sensitivity and 
high spectral resolution observations with the IRAM 30 m antenna, which allowed us 
to perform detailed analysis of the kinematics within one IRDC, G035. 39-00. 33. This 
cloud has been selected for its highly filamentary morphology and the presence of 
extended quiescent regions, characteristics of dynamical youth. We focus on the J = 
1 — > and J = 3 — > 2 transitions of N2H + , C 18 (1 — 0), and make comparison with 
SiO (2 — 1) observations and extinction mapping. Three interacting filaments of gas 
are found. We report large-scale velocity coherence throughout the cloud, evidenced 
through small velocity gradients and relatively narrow line widths. This suggests that 
the merging of these filaments is somewhat "gentle" , possibly regulated by magnetic 
fields. This merging of filaments may be r espon sible for the weak parsec-scale SiO 
emission detected by iJimenez-Serra et alj ([2010), via grain mantle vaporization. A 
systematic velocity shift between the N2H + (1 — 0) and C 18 (1 — 0) gas throughout 
the cloud of 0.18 ± 0.04 kms" 1 is also found, consistent with a scenario of collisions 
between filaments which is still ongoing. The N2H+ (1 — 0) is extended throughout the 
IRDC and it does not only trace dense cores, as found in nearby low-mass star-forming 
regions. The average H2 number density across the IRDC is ~ 5xl0 4 cm~ 3 , at least 
one order of magnitude larger than in nearby molecular clouds where low-mass stars 
are forming. A temperature gradient perpendicular to the filament is found. From 
our study, we conclude that G035. 39-00. 33 (clearly seen in the extinction map and in 
N2H + ) has been formed via the collision between two relatively quiescent filaments 
with average densities of ~ 5 x 10 3 cm~ 3 , moving with relative velocities of ~ 5 km s~ 1 . 
The accumulation of material at the merging points started > 1 Myr ago and it is still 
ongoing. 

Key words: stars: formation; ISM: individual objects: G035. 39-00. 33; ISM: 
molecules. 
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1 INTRODUCTION 

In order to disentangle the complex nature of massive (> 
8 Mq) star and stellar cluster formation, we must first under- 
stand the initial conditions under which these objects form. 
One way to do this is to, firstly, observe dense clouds that 
have the properties (masses, densities, temperatures etc.) 
required to produce massive dense cores (hereafter MDCs). 



2 Henshaw, Caselli, Fontani, Jimenez- Serra, Tan, Hernandez 




Figure 1. Average spectra of (Red) the isolated component of N 2 H+(1 - 0), (Black) N 2 H+(3 - 2), (Blue) C ls O (1 - 0), and (Green) 
SiO (2 - 1). The N 2 H+(1 - 0), (3 - 2) and SiO (2 - 1) intensities have been multiplied by factors of 2, 1.5 and 20, respectively. 



Secondly, it is preferentially required that the MDCs are lo- 
cated sufficiently far away from active star formation sites. 
This is to avoid complications due to, for instance, feedback 
effects (such as outflows, stellar winds or UV radiation), and 
to focus solely on 'pristine' MDCs. 

Within the past two decades, one particular group of 
objects has emerged as an excellent candidate to host the ini- 
tial conditions of star and stellar cluster formation, Infrared 
Dark Clouds (hereafter IRDCs). Discovered in the 1990's 
during mid-IR mapping of the Galaxy wit h the Infrared 
Space Observatory (ISO; IPerault et al. and the Mid- 

course Space Experiment (MffX; lEgan et al.1ll99ct ), IRDCs 
are high-extinction clouds, seen in silhouette against the 
diffuse mid-IR Galactic background. Studies have shown 
IRDCs are typically co l d (T < 25 K; e .g. IPillai et al.1 
l200d iPeretto et al.ll2pjltil , iRagan et al.l 120111 ). dense (n(H 2 ) 

10 5 M Q ), have large col- 
2 ^ 'e.g.lCarev et al.lll99Sl. 



> 10 5 cm~ 3 ), massive (~ 10 - 
umn densities (N(H 2 ) ^ 10 22 cm" 2 ) (e.. . 
Rathborne et al.1 l200rj . [§ imon et all l20od. Ivasvunina et al 



and have large d e uterium fractions (e.g . IPillai et al 



20071 . IChen et alJl201Ch . iFontani et ail (|201ll) found large 



N 2 D + /N 2 H + column density ratios towards a sample of high 
mass starless cores embedded in cold IRDCs (~ 0.3 — 0.7). 
IRDCs have a range of morphologies, and those having a 
filamentary structure are thought to be in an early stage of 
evolution, as describ ed by dynamical formation models of 
molecular clouds (e.g. | van Loo et afll2007l . iHennebelle et al.l 
l200Sl , lHeitsch et al.ll2009l ). 

The IRDC studied here, G035. 39-00. 33, has been se- 
lected to (i) have one of the most extr eme filamentary struc- 
tures in the iRathborne et al.l ([2006) sample of 38 IRDCs 



studied in millimetre continuum emission, (ii) have extended 
quiescent regions with relatively low star formation activity, 
and (iii) be relatively nearby (D = 2.9 kpc). We have fo- 
cussed specifically on the northern portion of the cloud. The 
southern portion of the cloud exhibits more tracers of star 
formation activity (8 /J,m and 24 fim point sources), and is 
thought to be in more of an advanced evolutionary stage. 
The portion over which we have mapped IRDC G035.39- 
00.33 conta ins three massive cores from the IRathborne et alj 
l|2006h and iButler fc Tar] |2009l . |2012| ; hereafter BT09 and 
BT12 studies). The most massive core (called H6 in BT09 
and MM7 in Rathborne et al. 2006) has a mass of 60 M 
within a radius of 0.15 pc (BT12) and there are no infrared 
sources at the mass surface density peak. Core H5/MM6 
hosts young stellar objects and has a mass of 36 Mq in- 
side a radius of about 0.13 pc, while H4/MM8 is starless, 
with a mass of 30 Mq inside 0.2 pc. iNguven Luong et alj 
(|201lh have also mapped IRDC G035.39-00.33 in the far- 
infrared with Herschel to investigate the star formation ac- 
tivity within the cloud. They detected a total of 28 cores, 
13 of which were considered to be MDCs (with masses > 
20 Mq, densities > 2 x 10 5 cm -3 ), potentially forming high 
mass stars. This study was carried out over the full length of 
the IRDC, and not just toward the northern portion selected 
here. 

Detailed mapping and analysis has recently been car- 
ried out towards the northern portion of IRDC G035.39- 
00.33. BT09 and BT12 used extinction mapping in the 8 /jm 
Spitzer-IRAC ba n d on a subsample of 10 IRDCs from the 
IRathborne et all (|2006l ) sample, providing a high angular 
resolution (~ 2") map of the mass surface density in all 
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Figure 2. Integrated intensity maps seen in contours, of (Left) 
N 2 H+ (1-0) isolated component and (Right) C ls O (1-0), over- 
laid on the mass surface density map smoothed to a 26" beam. 
Contours arc (Left) Wo- to 1.5K kms -1 , by 5a, where <r=0.05K 
kms -1 , (Right) 10ct to 4.5K kms -1 , by 5a, where cr=0.15K 
km s _1 . The dashed contour at 5.0 K km s _1 is added to show the 
location of the peak. The red and blue circles in the top left of each 
panel refer to the half-power beam widths of the N 2 H+ (1 — 0), 
and C ls O (1 — 0) maps (26" and 23", respectively). White crosses 
indicate locations of the cores identified by BT 12 following their 
detect ion in millimetre continuum emission by iRathborne et al] 
j2006h . 



cases. [jimenez-Serra et all (|2010l ) (hereafter, Paper I), pre- 
sented the first results from a dedicated molecular line study 
of this IRDC. They found widespread SiO (2 — 1) emis- 
sion covering a large area of the filament. This is sugges- 
tive of a large-scale shock caused by the recent collision 
of molecular filaments, which may have formed the IRDC. 
I Hernandez et all |201ll ) (hereafter, Paper II) measured sig- 
nificant levels of CO freeze-out acros s the I RDC, indicative 
of cold dense gas. iHernandez k, Tar] (|201ll 'l studied the dy- 
namics of the IRDC using 13 CO (1 — 0) from the Galactic 
Ring Survey (Jackson et al. 2006). Their results suggested 
the possible presence of dynamically important surface pres- 
sures around the cloud, perhaps indicating it was relatively 
young and had not yet settle d into virial equilibrium. How- 
ever, [Hernandez^£|al] (|2012h (hereafter, Paper III), revis- 
ited this question of virial equilibrium, now using higher 
resolution C ls O data and the improved extin ction map of 
BT12. In contrast to IHernandez fc Tanl (120111 ). they found 
the IRDC exhibited kinematics consistent with models of 
filamentary virial equilibrium, which requires the cloud to 
be older than one dynamical time, ~lMyr. 

In this paper (Paper IV of this series), we focus mainly 
on N 2 H+ (1 - 0), N 2 H+ (3 - 2) and C ls O (1 - 0). N 2 H+ 
is a known tracer of cold dense regions (> 10 4 cm -3 ), and 
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Figure 3. Integrated intensity contours of N 2 H+ (3 — 2), overlaid 
on the mass surface density map smoothed to a 9" beam. Con- 
tours are 3<r to 5.52Kkms — 1 , by la, where cr=0.46K kms - 1 . 
The dotted contour is the 2a contour. White crosses indicate 
the positions of the massive cores (BT12). Colored symbols are 
4.5 fim (green squares), 8 fim (red circles) and 24 fim (red trian- 
gles) sources found in the direction of the IRDC (see Figure 1 of 
Paper I for more details on these sources). 



in low-mass star-for ming regions it there fore typically traces 
the dense cores (e.g. lCaselli et al"1 [2002a). C ls O is a tracer of 
more diffuse material (^ 10 3 cm -3 ). Due to the fact that it 
is readily abundant in molecular clouds, C 18 gives more of 
an indication of the general morphology of the whole cloud. 
Unlike 12 CO and 13 CO, the relatively less abundant C 18 
is typically optically thin in molecular clouds, making it an 
ideal diagnostic tool to study the kinematics of the lower 
density material. In addition to these tracers we also include 
the mass surface density derived by BT12 as a tracer of the 
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dense material in the filament. The SiO (2 — 1) emission, 
studied extensively in Paper I, is also used to compare our 
kinematic results, spatially, with the shocked gas within the 
filament. Finally, the J = 3 — > 2 transition of N2H + is used 
together with the J = 1 — » transition to estimate the 
number density across the filament. 

Details of the observations can be found in Section[2l 
with the results (integrated intensity maps, channel maps, 
kinematics, column and number density) in Section[3] A dis- 
cussion of our results can be found in Section[4] and con- 
clusions are in Section[5] AppendixlAl outlines the fitting 
method used to extract information from this data, whereas 
Appendix B reports tables with the Gaussian and hyperfine 
structure fits in selected positions across the filament. 



2 OBSERVATIONS 

The N 2 H+ (1 - 0), N 2 H+ (3 - 2), C ls O (1 - 0), and SiO 
(2 — 1) observations were mapped with the Instituto de Ra- 
dioastronomi'a Milimetrica (IRAM) 30-m telescope at Pico 
Veleta, Spain. The N2H + observations were carried out in 
August 2008 and February 2009, C ls O (1 - 0) in August 
2008, and SiO (2 — 1) observations in December 2008 and 
February 2009. The large-scale images were obtained with 
the On-The-Fly (OTF) mapping mode using the offsets 
(1830", 658") for the 2008 observations, and (300", 0") for 
observations carried out in 2009 (i.e. for part of the N2H 4 " 
and SiO data, whose emission is concentrated in the IRDC), 
as off-positions. The central coordinates of the maps are 
a(J2000) = 18' i 57 m 08 s , 5(J2000) = 02°10'30" (I = 35°. 517, 
b = -0°.274). All observations were carried out using the 
old ABCD receivers. The VErsatile SPectrometer Assem- 
bly (VESPA) provided spectral resolutions between 20 and 
80 kHz. Information on the beam sizes, frequencies and ve- 
locity resolutions can all be found in Table[T] Typical sys- 
tem temperatures range from 100-200 K for all line obser- 
vations except for N2H + (3 — 2) which has a value of Tsys 
~ 750 K. All intensities were calibrated in units of antenna 
temperature, T^_. To convert these intensities into units of 
main-beam brightness temperature (assuming a unity filling 
factor), Tmb, we have used the beam and forward efficien- 
cies listed in TableQ] Saturn was used to calculate the fo- 
cus, and pointing was checked every 2 hours on G34.3+0.2. 
The data were calibrat ed with the chopper-wheel technique 
(jKutner fc Ulichl Il98lh . with a calibration uncertainty of 
~20%. 

As in Paper III, we utilize the 8 fj,m extinction de- 
rived mass surface densit y map of BT12, as modified by 
iKainulainen fe Tanl l|2012l ) (hereafter, KT12) to include cor- 
rections for the presence of the NIR extinction derived IRDC 
envelope. The mass surface density map has 2" resolution, 
which we have smoothed to 26" (unless otherwise stated), so 
as the data is of comparable resolution to the N2H + (1 — 0) 
and C ls O (1 — 0) observations. 



3 RESULTS 

3.1 Spectra and integrated intensities 

The spectra of the isolated N 2 H + (1 — 0) hyperfine compo- 
nent (Fi, F = 0,1 -> 1,2), the N 2 H+ (3-2), the C 18 (1-0), 



and the SiO (2 — 1), averaged across the whole mapped area, 
are shown in Figure[T] The spectra are plotted in units of 
main beam brightness temperature, with the N2H + (1 — 0), 
N 2 H+ (3 - 2) and SiO (2 - 1) being multipled by factors of 
2, 1.5 and 20, respectively, for comparison. The four molec- 
ular gas tracers have different centroid velocities, with the 
highest density tracer, N2H + (3 — 2), peaking at the largest 
velocity, followed by N 2 H+ (1 - 0), C ls O (1 - 0) and SiO 
(2 - 1). The SiO (2 - 1) shows the broadest width, with 
emission present up to velocities of about 50kms _1 . The 
C ls O line encloses the emission of the other species (ex- 
cept for the high velocity emission of the SiO (2 — 1)) and 
presents an asymmetric profile indicative of multiple velocity 
components along the line of sight. Although shifts in veloc- 
ities relative to the SiO line are expe cted given that SiO 
typic ally traces shocked material (e.g. Ijimenez-Serra et al.1 
I2009T ). the significant shifts in centroid velocities observed in 
the more quiescent molecular cloud tracers C ls O and N2H + 
repres ent a clear difference from low-mass star -forming re- 
gions l|Walsh et alj|2004 lHacar fc Tafallall201ll ). 

Figure[2] presents integrated intensity maps of the J — 
1 — > transitions of both N2H + (Red contours; Fi, F = 0,1 
— > 1,2 hyperfine component), and C ls O (Blue contours), 
convolved with the respective beams (see Table[T} and using 
a pixel size of half the respective beam sizes. In each map, 
the spectra have been integrated between 40 and 50 kms , 
to cover all the emission (see Figure[T]|. The N2H + (1 — 0) 
and C 18 (1 - 0) maps have been superimposed on the 
mass surface density map smoothed to 26" angular reso- 
lution. Overlaid are white crosses to indicate the location 
of the thr ee massive cores seen i n the millimetre continuum 
emission (Rat hborne et a l. 2006) and in the unsmoothed ex- 
tinction map (|Butler fc Tanll2009 ) . The integrated intensity 
map of the N2H + (3 — 2) emission is shown in Figure(3] over- 
laid on the mass surface density map smoothed to the IRAM 
9" beam at the frequency of this transition. 

From Figures[2] and |3] we find that the N 2 H+ (1 - 0) 
and (3 — 2) emitting regions are extended across a large por- 
tion of the filament. This is different from low-mass star- 
forming regions, where N2H + mainly traces the dense core s 
ijCaselli et al.ll2002al . lAndre et al1l2007l . iFriesen et aDl2010T ) , 
whereas filaments can be seen with lower density tracers 
such as C ls O (1 - 0) (|Mizuno et al.lll995l . lHacar fc Tafallal 
Hon]). Given that the critical density of the N 2 H+ (1-0) 
transition is ~10 5 cm~ 3 , this is a clear indication that the 
overall number density of this IRDC is similar to the aver- 
age density of nearby low-mass dense cores (i.e. at least one 
order of magnitude higher than the average density of the 
clouds within which nearby low-mass cores are embedded, 
<10 3 cm~ 3 ; e.g. Pineda et al. 2008). N 2 H+ (3 - 2) peaks 
towards two out of the three massive cores mentioned above 
(white crosses in Figure|3J). The peak in both transitions of 
N 2 H + is seen to be coincident with the position of core H6, 
whereas the C ls O (1 — 0) peaks away from this (within ~ 1 
beam size). This is due to the freeze-out of CO, most promi- 
nent towards core H6 (Paper II). Emission from both N2H + 
transitions traces the morphology of the extinction fairly 
well, whereas the emission from the C ls O (1 - 0) covers a 
larger area. 

Figure^ shows the total integrated intensities of the 
J = 1 -> (Fi, F = 0,1 — > 1,2) and J = 3 -> 2 transi- 
tions of N 2 H+, C 18 (1 - 0), and SiO (2 - 1) as a function 
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Figure 4. Integrated intensity as a function of mass surface density for the (Red) N2H + (1 — 0) isolated component, (Blue) C ls O (1 — 0), 
(Green) SiO (2 — 1) (all smootcd to 26"), and the 9" resolution (Black) N2H+ (3 — 2). The N2H+ integrated intensity spectra within 
the starless core H6 are displayed in Cyan. 



of mass surface density. Intensities are integrated over a ve- 
locity range of 40-50 kms -1 . To highlight the properties of 
the most massive core in the region, N2H + data from core 
H6 have different color f]. Significant correlation (regression 
value, r-value = 0.8) is found between N2H + (1 — 0) inte- 
grated intensity and mass surface density, over the full cloud 
(including core H6). This is to be expected, as the N2H + is 
tracing the dense gas. A significant correlation is also found 
across core H6 (r-value = 0.6), despite the fewer data points. 
For the N2H + (3 — 2) we have compared the line integrated 
intensity and mass surface density, both smoothed to 9". 
The correlation for N2H + (3 — 2) is however not as strong 
(r-value = 0.5), over the full cloud, probably due to the 
fact that this line is not just sen sitive to number d ensity 
but also to temperature (see e.g. Pa gani et al J 12007), and 
warm gas may not be correlated with the mass surface and 
number density (see Section [33]) . A slightly stronger corre- 

1 In order to separate this core from the rest of the data we 
have taken the 5cr level of the N2H+ (1 — 0) and (3 — 2) intensity 
(between -20.0 s: Act < 20.0 and -10.0 < A<5 ^ 50.0), and attribute 
this to the core. 



lation is present toward core H6 (r-value = 0.55), where the 
line is brightest. There is also a strong correlation (r-value 
= 0.8) between the C ls O (1 — 0) integrated intensity and 
the mass surface density. A weaker correlation (r-value = 
0.45) between the SiO (2 — 1) integrated intensity and the 
mass surface density is observed. This, along with the fact 
that the centroid velocities of the SiO and N2H + transitions 
are not coincident (see Figure]!}, suggests that the narrow 
SiO (2 — 1) emission is mainly associated with the lower 
density gas traced by C ls O (uh ~ 10 3 cm~ 3 ), despite the 
relatively high critical density of the transition (n cr ~ 10 5 
cm -3 ). Therefore, we expect the narrow SiO (2 — 1) lines 
to be subthermally excited, which would explain their low 
brightness (Paper I). 



3.2 Kinematics 

3.2.1 Multiple velocity components 

Figure[S] displays the individual N 2 H+ (1 - 0) (Fi, F = 0,1 
— > 1,2) and C ls O (1 — 0) spectra over the full extent of 
the region mapped in N2H 4 " (1 — 0). This is overlaid on 
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Figure 5. Spectra of the (Red) N 2 H+ (1 - 0) (Fi, F = 0,1 -> 1,2) and (Blue) C ls O (1 - 0) at all positions in the cloud, overlaid on the 
mass surface density map of KT12. Spectra are shown between 41-49 km s _1 and -0.1-1.2 K. The intensity of the N2H+ (1 — 0) has been 
multiplied by a factor of 2 for clarity. The C ls O (1 — 0) has been smoothed to 26", and regridded this to 13", to make direct comparison 
with the N 2 H+ (1-0) data. 
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Table 1. Observed molecular transitions, line frequencies, spectral resolutions, telescope beam sizes, beam correction factors for the 
ABCD receivers at the IRAM 30 m telescope. 



Transition 


Frequency 


Spectral 


Beam Size 


Beam 


Forward 




(MHz) 


Resolution (kms -1 ) 


(arcsec) 


Efficiency 


Efficiency 


SiO J = 2 -> 1 


86846.9600° 


2.68X10- 1 


28 


0.77 


0.98 


N 2 H+ J = 1^0 


93173.7637 b 


6.28xl0~ 2 


26 


0.76 


0.98 


C 18 OJ = 1^0 


109782. 1730 c 


5.33xl0~ 2 


23 


0.73 


0.97 


N 2 H+ J = 3^2 


279511.8320 6 


4.19xl0~ 2 


9 


0.49 


0.90 



Pagani et alj d2009 



Cazzoli et al] J2003 




60 0-60 
Aa (arcsec) 



Figure 6. G035. 39-00.33 channel maps of (Red) N 2 H+ (1 - 0) (Fi, F = 0,1 -> 1,2) and (Blue) C ls O (1 - 0), overlaid on the mass 
surface density map. Intensity maps integrated from 43 kms -1 to 48kms _1 and in steps of 0.2kms _1 . Contours for N2H+ (1 — 0) in 
all cases are from 3cr, increasing in steps of 3<r (0.02 Kkms -1 ). Contours for C 18 (1 — 0) in all cases are from 3<r, increasing in steps 
of 3cr (0.04 Kkms" 1 ). 



top of the mass surface density map, smoothed to 26", to 
match the angular resolution. As previously noted in Paper 
I, multiple velocity components appear at various positions 
within the cloud (e.g. around core H6 and toward the South). 
This is more apparent in C ls O (1 — 0) as it traces the more 
abundant, less dense gas along the line of sight. However, the 
line morphologies of the C ls O (1 — 0) are echoed in the South 
and Central positions in N2H + (1 — 0). The fact that the 
profiles of the N2H + (1 — 0) line are similar to the optically 
thin C ls O (1 — 0), means that we are indeed seeing multiple 
velocity components along the line of sight and not optical 
depth effects. In particular, the N2H + (1 — 0) spectra around 
the peak of core H6 (see offsets (-7.5, 8), (-7.5, 20), (5.5, 8), 



and (5.5, 21)), with a stronger blue pe ak typically observe d 
in collapsing low-mass dense cores (e.a ACrapsi et aZ]|g00.Sj) . 
are just the superposition of different velocity components. 

A more detailed view of the gas distribution can be 
found in Figure[B] Here we present the channel maps of the 
emission seen in N 2 H + (1 - 0) (Fi, F = 0,1 -> 1,2) and 
C ls O (1 — 0), between 42-48 kms -1 , in velocity increments 
of 0.2 kms -1 , superimposed onto the mass surface density 
map smoothed to 26". From the figure, it is evident that the 
additional components seen in the spectra of Figure[5] have 
varying morphologies (see also Paper I for a similar analysis 
using C 18 (2 — 1) data). The first of these components (Fil- 
ament 1) can be seen between 42-44 kms -1 . In C 18 Q (1-0) 
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Figure 7. Position- Velocity diagrams of both (Red) N2H+ (1 — 0) (Fi, F = 0,1 — > 1,2) isolated component and (Grey-scale) C ls O 
(1 — 0). Cuts are taken from North to South through the cloud at each offset right ascension (which can be found in the top left of each 
panel). Emission in N 2 H+ (1 - 0) and C ls O (1 - 0) is plotted from 3<r (<r = 0.14, and for N 2 H+ (1 - 0), a = 0.06). Contours increase in 
3<r steps. Vertical dotted lines refer to the centroid velocities derived from the GGF (see Appendix fXJ . Horizontal dashed lines at offsets 
Aq = -7.5 and 5.5, refer to the positions of core H6 (2, 20), and core H5 (6, -86). 



the component is seen to run from North-East to South- 
East, curving to the West toward the centre of the map. 
The N2H + (1 — 0) emission in this component is restricted to 
the most southern portion of the cloud. The gas distribution 
changes shape between 44-46 km s (Filament 2), peaking 
at the positions of the massive cores. Here, the N2H + (1 — 0) 
follows the morphology of the C 18 (1 — 0) more closely, and 
is extended over the whole filament, indicating the presence 
of widespread dense gas. Both species follow the shape of 
the extinction map, from North- West to the South, curving 
slightly to the East. Finally, another component is seen be- 
tween 46-48 kms -1 (Filament 3). This emission is similar in 
morphology to the previous component. However, the emis- 
sion is restricted to the northern portion of the cloud, and 
is not seen in the South. Given that this component has a 
well defined feature in both the spectra (Figure[5]), and the 
channel maps (Figure^ , we consider this to be a separate 
filament. 

Another way to inspect the gas kinematics across the 



cloud is through position-velocity (PV) diagrams. Figure[7] 
shows the PV diagrams of both species (N 2 H + (1 — 0) (Fi, F 
= 0,1 — ¥ 1,2) in red, and C ls O (1 — 0) in grey scale), slicing 
the cloud from North to South, at each offset of right as- 
cension in our N2H + (1 — 0) map. Vertical dotted lines refer 
to the velocities of the individual components seen in both 
species as deduced by the Guided Gaussian Fit (GGF) of 
the average spectrum (see Appendix A for further details). 
At the most westerly point of the cloud (Aq = -33.5"), the 
N2H + (1 — 0) emission follows the C 18 (1 - 0) emission 
in both the North and South, except at AS ~ 30", where 
we can see an additional component in C 18 (1 — 0) which 
is part of Filament 1. There is no significant velocity struc- 
ture, except for a velocity change of about 1 kms -1 between 
the northern and southern emission on a scale of ~ 1 pc, 
which corresponds to a very small velocity gradient of ~ 
1 kms" 1 pc" 1 . Such a small g radient is comparab le to those 
found in low-mass cores (e.g. ICaselli et al.1l2002al ). but here 
it is seen over a significantly larger extent. 
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Figure 8. Vlsr ar >d Av maps of Filament 2, as derived from the Guided Gaussian Fits (see Appendix [Aj . Maps on the left represent 
the V LS r of (Left) N 2 H+ (1 - 0) and (Right) C ls O (1 - 0). Maps on the right represent the linewidth, Av of (Left) N 2 H+ (1 - 0) and 
(Right) C 18 (1 — 0). Black crosses indicate the positions of the massive cores from BT12. 



As we move towards the centre of the cloud, the emis- 
sion becomes stronger (max intensity of N2H + (1 — 0) (Fi, 
F = 0,1 -» 1,2) at -33.5" = 0.82 K km s" 1 compared to 
1.25 K kms" 1 at 5.5") and spreads over a larger velocity 
range. This can be seen towards offsets Aa = -20.5" and 
-7.5", where an additional component at higher velocities 
(Filament 3) appears (most prominent between 0" < AS 
< 50"). Core H6 lies between -7.5" < Aa < 5.5" at AS = 
20". Here all three filaments are seen, with Filament 2 being 
the most prominent in N2H + (1 — 0). The channel maps of 
Figure[6] also show that these features coincide spatially at 
the position of core H6, suggestive of filament merging (see 
below) . 

3.2.2 Centroid velocity and line width 

Figure[5] shows the centroid velocity map of Filament 2 in 
N 2 H+ (1 - 0) (Fi, F = 0,1 -> 1,2) and C 18 (1 - 00. 
These velocities are derived from the GGF (Guided Gaus- 
sian Fit) method (Appendix A). The velocity field appears 
coherent, except for a velocity change of ~ 0.5 kms" 1 at the 

2 Both the low and high velocity components (Filament 1 and 3, 
respectively) do not cover a large enough area in N 2 H + (1 — 0) 
to see any significant velocity gradient, so we do not consider the 
other two filaments in this analysis. 



northern edge of core H6. The change in velocity from ~ 
45.9 kms" 1 (North) to ~ 45.4 kms" 1 (South) occurs within 
~ one beam width, corresponding to a local velocity gra- 
dient of ~lkms~ 1 pc~ 1 . Same conclusions are drawn from 
the C 18 (1 — 0) line, although the velocity field appears 
more complex south of core H6. Thus, the abrupt (but rela- 
tively small) velocity change is present in the dense material 
as well as in the lower density envelope, suggesting that it 
originates from the larger scale. 

Figured] also shows the N 2 H+ (1 - 0) (Fi, F = 0,1 -> 
1,2) and C ls O (1 - 0) line width maps of Filament 2. The 
first thing to note is that line widths are dominated by non- 
thermal motions. In fact, adopting an a verage kinetic tem- 
perature of 15 K (e.g. iRagan et al.ll201ll ). the thermal width 
for both species is 0.15 kms" 1 . Thus, the observed widths 
are between 3 and 15 times the thermal width, with the nar- 
rowest widths found toward the North. On average, C 18 
(1 — 0) and N 2 H + (1 — 0) line widths have similar values. 
This is in striking contrast to low-mass star-forming regions, 
where lines o/C ls O (1—0) typically show broader line widths 
compared to those measure d in high density tracer s (by a fac- 
tor between 1.5 and 2; e.a. \Fuller & Mveri\l99^{) . However, 
this is consistent with our finding that N 2 H + is widespread 
across the filament and not preferentially tracing dense cores 
as in low-mass star- forming regions (see Figure 2). In both 
species, the line widths get narrower towards the northern 
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portion of the cloud and (to a lesser extent) toward the star- 
less core H6, indicating the presence of relatively quiescent 
(and probably pristine) gas. 

Figure^ displays the V LS r of N 2 H+ (1 - 0) (Fi, F = 
0,1 -> 1,2) and C 18 (1 - 0) as a function of mass surface 
density. The three filaments are clearly seen as well sepa- 
rated velocity components (as determined by the Gaussian 
fits). The average velocities of the individual components 
differ slightly between the two species. In Filament 1, we 
find a mean V LSR of 42.95 (0.17) kms" 1 in N 2 H+ (1 - 0), 
and 43.65 (0.12) kms -1 in C ls O (1 - 0). Filament 2 has a 
mean velocity of 45.63 (0.03) kms" 1 in N 2 H+ (1 - 0), and 
45.40 (0.03) kms" 1 in C 18 (1 - 0). Filament 3 has a mean 
velocity of 46.77 (0.06) kms -1 in N 2 H+ (1 - 0), and 46.76 
(0.05) kms" 1 in C ls O (1 - 0). We note that at mass sur- 
face densites below 0.06 gem" 2 , the C ls O filaments span 
a broader range of velocities. Moreover, the relative veloc- 
ity between Filament 2 and 1 is larger than that between 
Filament 2 and 3. This could be at the origin of the ob- 
served velocity shift between N 2 H + and C 18 in Filament 2 



(see Figure[5]and sect ions 1 3 . 2 . 3 1 and l4l for more details). The 
abrupt change in velocity observed in Filament 2 (Figure[7] 
and [8]) is seen here as the "gap" between the two groups of 
blue points. It is interesting to see that the velocities asso- 
ciated with the massive starless core H6 are filling this gap, 
suggesting that this core has probably formed at the interface 
of material moving at different velocities. 



3.2.3 The N 2 H+- C 18 velocity shift 

Figures[l0] and [TT] make clear that the velocity shift seen 
between N 2 H+ (1-0) and C ls O (1-0) is the result of large- 
scale kinematics. Figure[l0] is a map of the N 2 H+ - C 18 
velocity shift measured across Filament 2, with the contours 
of SiO (2—1) overlapped on top. The N 2 H + (1—0) emission is 
largely red-shifted with respect to the C 18 (1 — 0) emission, 
and we report a mean velocity shift of 0.18 ± 0.04 kms" 1 , 
across the whole filament. The largest velocity shifts are seen 
to the North and South of the cloud. North of offset AS = 
86.0, the average velocity shift is 0.22 (0.04) kms" 1 . South 
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of offset AS = -70.0", the velocity shift is 0.26 (0.04) kms -1 . 
Between these offsets, the velocity shift is 0.13 (0.04) kms -1 . 
No correlation is found between the velocity shift and the 
SiO (2—1) integrated intensity. The velocity shift is also clear 
in Figure fTTI which shows the velocity difference between 
C ls O and N2H + as a function of mass surface density. Also 
in this case, no correlation is found, but the shift is clearly 
present all across the filament. 

The velocity shift is not constant along the filament. 
Figure [T2l shows the N2H + - C ls O velocity difference as a 
function of declination along the four central strips at fixed 
right ascension offsets -7.5", 5.5", 18.5" and 31.5". From 
this figure it is clear that the shifts in the North and South 
have similar magnitudes, while variations are present in be- 
tween. The variations along the four strips appear similar 
up to the H6 core, where the outer strip at Aa = -7.5" 
(and to a lesser extent the strip passing through the peak of 
core H6) show different behavior. This is the region where 
the abrupt change in velocity has been found (see Figure|S| . 
Those regions where the dense gas and its envelope appear 
to be moving with the same velocity are present at AS = 
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Figure 11. Velocity shift between N 2 H+ (1-0) and C ls O (1-0) 
as a function of mass surface density. The right of the plot shows 
a histogram of this velocity shift information. Errors in the mass 
surface density represent an assumed uncertainty of 0.007 gcm~ 2 . 



-40" in the two western strips and at the peak of core H6. 
However these points of minimum shift are surrounded by 
sharp variations, especially evident toward core H6. Discus- 
sion on a possible origin of this velocity shift can be found 
in section 13] 

3.3 The C ls O/N 2 H + integrated intensity ratio 

Tafall a~fc Santiago! (|2004l ) used CO depletion as a chemi- 
cal clock, to indicate an evolutionary sequence of individual 
cores. This is achieved by taking two species that show dif- 
ferent amounts of depletion at high densities. C ls O freeze- 
out becomes heavy at densities upwards of a few 10 4 cm -3 , 
whereas N 2 H+ (1-0) remains in the gas phase up to dens i- 
ties > 10 5 cm" 3 (e.g.. lCaselli et al.lll999l . lTafalla et afl2002D . 
Moreover, N2H + takes longer than C ls O to form, due to 
the slower neutral-neutral reactions involv ed in the produc- 
tion o f N 2 , the parent species of N2H + (e.g. iHilv-Blant et alj 
2010). Therefore, an early stage of evolution would be char- 
acterised as having an undepleted C 18 and low abun- 
dances of N2H + . As the cores follow their evolution toward 
star formation, their density increases and C ls O starts to 
frcczc-out. Thus, with time, the CO abundance falls while 
N2H + abundances reach equilibrium and maintain relatively 
large values up to volume densities of ^ 10 5 cm -3 , as the 
freeze-out of CO initially boosts the production rate of 
N2H + (given that CO is an important destruction part- 
ner of N2H + ). In light of this, we can use the method of 



iTafalla fc Santiago! (120041 ) and take the ratio of the inte- 
grated intensities of C ls O (1 - 0) and N 2 H+ (1 - 0): 



R : 



7(C ls O) 



J(N 2 H+) 



(1) 



to define a chemical clock. Chemically young gas is charac- 
terized by values of R > 1. Figure [T3l displays the map of 
the C ls O/N2H + intensity ratio integrated between 44 and 
46 kms -1 , where Filament 2 is more prominent. To obtain 
this map, we first smoothed the C ls O (1 - 0) to 26", the 
angular resolution of the N2H + (1 — 0) map, and regridded 
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the C 18 (1 - 0) map. The C 18 0/N 2 H+ integrated inten- 
sity ratio has a value < 1 throughout the whole filament, 
except for the northern regions. The bulk of the filament 
shows large volume of N2H + with respect to C 18 0, and 
this suggests that the gas is in a chemically evolved state. 
This is clearly seen in Figure fTH where R is plotted ver- 
sus the mass surface density: no R values larger than 1.0 
are found at £ ^ 0.09 gem -2 . Given that the chemical time 
scale is inversely dependent on the number density, the ob- 
served advanced chemical evolution is probably just the result 
of the relatively large average density of the filament when 
compared to low-mass star-forming regions (see Section l3~5")) . 
Only in the North, the gas appears relatively less chemically 
evolved, suggesting the presence of lower density material. 



3.4 Column Density 

The N2H + and C ls O column densities have been calculated 
from the integrated intensity of the N2H + (1 — 0) Fi, F = 
0,1 — s> 1,2 hyperfine component (divided by its statistical 
weight) and C 18 (1 — 0) emission. The assumption here 
is that both lines are optically thin. This has been verified 
for the N2H" 1 " isolated compon ent (with statistical weight 
of ~ 0.11111; ICaselli et aflll995h with the HFS fitting tech- 
nique (see Appendix B) and for C ls O (1 - 0), Paper II). To 
measure the column densities, we then used equation (A4) 
in Caselli et al. (2002). 

The average N 2 H + column density across the IRDC is 
N(N 2 H+) = (1.14 ± 0.11) x 10 13 cm -2 , peaking in the 
South at offset (5.5", -122"), with a value of (2.33 ± 0.25) 
x 10 cm - . The minimum value of N(N 2 H + ) is located 
in the North at offset (5.5", 125"), where N(N 2 H+) = 
(2.33 ± 0.69) x 10 12 cm -2 . These column den sities are sim- 
ilar to those measured toward low-mass cores ( Case lli et alj 
l2002al ). indicating that the average physical and chemical 
conditions of the whole filament traced by N2H + (1 — 0) are 
similar to those of the well-known dense cores in nearby low- 
mass star-forming regi ons, in agreement with the study of 
IVasvunina et alj < |201lh . 

The average C ls O column density throughout the fila- 
ment is N(C ls O) = (3.82 ± 0.04) x 10 15 cm -2 , peaking at 
position (-10.0, 34.0), North- West of the extinction peak of 
core H6 with a value of (5.95 ± 0.72) x 10 15 cm" 2 . This is 
consistent with Paper II, that found evidence of CO freeze- 
out toward the center of H6. 

We find a significant correlation (r-value = 0.8) between 
the column density of N2H + (1 — 0) and the mass sur- 
face density (as expected given the strong correlation found 
in Figure[4] between the N 2 H + (1 — 0) integrated intensity 
and the mass surface density): N(N 2 H + ) = (1.5±0.2) x 10 12 
+ (9.8±0.3) x 10 13 E. This correlation implies a constant 
fractional abundance of N 2 H + (w.r.t. H2 molecules) of 
~4xl0 -10 , simi lar to the value found in low-mass cor es (a 
few times lQ- 10 : ICaselli et al.ll2002bl . lCrapsi et"ai1l2005T l and 
in the quiescent mass ive starless core in Orion (Orion B9; 
iMiettinen et al.ll2012l ). but lower than the value measured 
in the direction of m assive IR-dark clumps (a few xl0~ 9 ; 
IMiettinen et al.ll201lf ). Quiescent dense core conditions once 
again appear to be appropriate in describing G035. 39-00. 33. 



3.5 Number density and kinetic temperature 

Observations of the J = 1 — > and J = 3 — > 2 rotational 
transitions of N 2 H + can constrain the number density, if the 
gas temperature is known (see Section [3~l"T) . No gas temper- 
ature measurements have been made yet across G035.39- 
00.33, so we assume a constant temperature of 15 K across 
the cloud, b ased on the kinetic temperature deriv ed in other 
IRDCs fe.g. lPillai et al.ll2006l , feagan et al.ll201lf) . and close 
to the dust tempera t ure m easured across the filament by 
iNguven Luong et alj (|201 lh using Herschel data. 

To calculate the nu mber density, the non-LT E radiative 
transfer code RADEX (|van der Tak et al.l [20071 was used. 
The N 2 H + (3 — 2) map was smoothed to 26" resolution, to 
match that of the N2H + (1 — 0) data. Given the complexity 
of the spectra and the fact that the (3 — 2) transition does 
not have isolated hyperfine components, we fitted a single 
Gaussian peak to the isolated hyperfine component of the 
N2H + (1 — 0) irrespective of the fact that we may see mul- 
tiple components along the line of sight. In addition to this, 
we performed a HFS fit assuming optically thin conditions 
to the N2H + (3 — 2) lines in CLASS, which gives, as one of 
its outputs, the linewidth. The average linewidth between 
the N2H + (1 — 0) and (3 — 2) profiles was chosen as input 
parameter for the RADEX calculation. The column density 
input to RADEX was calculated using the integrated inten- 
sity of N2H + (1 — 0) over the velocity range 40-50 km s _1 
(first using the isolated component and then scaling by the 
statistical weight). As a final assumption we then used a 
value of 15 K for the kinetic temperature of the gas. The 
number density was then varied between lxlO 3 cm -3 and 
1x10 s cm -3 and a minimisation technique was used to find 
the best fit to the observed brightness temperature ratio 
of the J = 3 — > 2 and J = 1 — ¥ lines. Figure [T5l shows 
the result of this minimization technique, i.e. the number 
density map of G035. 39-00. 33. The average H2 number den- 
sity across the filament (measured with a 26" beam, corre- 
sponding to a linear size of 0.36 pc) is ~5xl0 4 cm -3 , with 
peaks close to the three massive cores and East of core H6. 
This value is roughly a factor of 10 larger than that found 
using extinction mapping in Paper III. Possible causes of 
this difference are unresolved clumping and the fact that 
the analysis done with the N2H + data only includes those 
points where the 3-2 line is observed, therefore it is biased 
toward higher densities (as the 3-2 line has a critical density 
of 3xl0 6 cm" 3 ). 

To check possible temperature variations across the fil- 
ament, we fix the number density to the average value found 
with the previous method and find the kinetic temperature 
map. Figurc[T6l shows the kinetic temperature map obtained 
from the N 2 H+ (3-2)/(l-0) intensity ratio and the RADEX 
code, assuming a constant number density across the fila- 
ment of 5xl0 4 cm -3 (the average value found in Figure [15)) . 
A temperature gradient is present perpendicular to the fil- 
ament, with the largest values (^40 K), close to the central 
SiO peak (see Figure [10}. Although this analysis should be 
taken with caution, given the simplistic assumption of con- 
stant density, we speculate that the observed temperature 
increase may be linked to the filament merging and the con- 
sequent shock which has produced the extended SiO emis- 
sion (Paper I). In fact, the extinction in these regions is still 
large enough to exclude any temperature enhancement due 
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Figure 12. Velocity shift between N2H + (1 — 0) and C 18 (1 — 0) as a function of offset in declination along four strips of constant 
right ascension (see legend in top right). The blue line refers to the strip passing closer to the center of core H6. Horizontal dotted lines 
correspond to the mean values of the velocity shift, 0.18±0.04 km s — 1 , and the zero line. Vertical dotted lines indicate the positions in 
offset declination of the three massive cores from BT12. 



to the interstellar radiation field. Independent temperature 
measurements (e.g. from the inversion transitions of ammo- 
nia) and higher angular resolution observations are needed 
to put more stringent constraints on the temperature and 
density structure across G035. 39-00. 33. 



4 DISCUSSION 

4.1 What is causing the N2H 4 
shift? 



C 18 velocity 



The velocity shift between the N 2 H+ (1-0), and C 18 (1-0) 
reveals relative motions between the dense gas (traced by 
N2H + ) and the surrounding less dense envelope (traced by 
C 18 0). This shift has been sear ched for but not found in low- 
mass star- forming regio ns (e.g .lWalsh eTa l. 2004, Kirk et all 
120071 . lHacar fc Tafallall201ll ). lHacar fc Tafallal (|201ll ) stud- 
ied the centroid velocities of N 2 H+ (1 - 0) and C ls O (1-0) 
in L1517, and found a good match between the two tracers. 
They use this to say that there are no significant motions be- 
tween the different density regimes of the cloud, namely the 



dense cores as traced by N2H + (1 — 0) and the surrounding 
envelope traced by C 18 (1 — 0). They go on to conclude 
that this implies that the velocity structure of the core is 
therefore not an intrinsic property of the core, but a result 
of large-scale motions on filamentary scales. In the case of 
G035.39-00.33, N 2 H+ (1 - 0) traces the dense cores as well 
as the filament over a distance of ~ 4 pc, similar to the C 18 
(1 — 0) line, so that the observed velocity shift represents a 
puzzle which needs to be investigated in more detail. 

The channel maps in Figure|B] clearly show that signifi- 
cant N 2 H + (1 — 0) emission is present in a narrower range of 
velocities compared to the C ls O (1 — 0) emission. In partic- 
ular, at the lowest velocities identified with Filament 1 (be- 
tween 42 and 44kms _1 ), N2H + is practically undetected 
along most of the filament, whereas Filament 3 (between 
about 46 and 48kms _1 ) is also seen in N 2 H + (1 — 0) up to 
47kms~ 1 . If the three identified filaments are interacting, 
the velocity field within the denser Filament 2 will then be 
affected by Filament 1 and 3 in the case of C ls O (1 — 0), but 
only by Filament 3 in the case of N2H + (1 — 0). This would 
imply a blueshift of the C ls O line, as we indeed observe 
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(see Figure fTTj) . The filament interaction is also evident in 
Figure[9] where the C ls O panel shows less distinct velocity 
components, and Filament 2 is split in two velocity compo- 
nents parallel to Filaments 1 and 3, with the densest region 
in the IRDC (core H6) in between. The abrupt change in 
velocity seen at the edge of core H6 (Figure|SJ) could also be 
produced by filament merging, with the consequent forma- 
tion of a dense core at the intersection. Finally, Figure [T2l 
displays a variation of the observed velocity shift between 
the C 18 (1 - 0) and N 2 H+ (1 - 0) lines along Filament 2, 
which can be explained by inhomogeneities along the inter- 
acting Filaments 1 and 3. Indeed, looking at Figures^ [TO] 
and[JJ3one can see that the maxima in the velocity shift seen 

in Filament 2 correspond to local maxima oftheC 18 (1-0) | Caselli et all 19971. 
intensity in Filament 1. At these positions the densities are 
expected to be higher (further indicated by the presence of 
N2H + (1 — 0)) and thus the interaction more prominent. 

The above observational evidence suggests that Fila- 
ment 2 may be the result of the merging of Filaments 1 and 
3, as sketched in Figure[l7] This process of filament merging 
is still ongoing, with material from Filament 1 still mov- 
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Figure 14. Integrated intensity ratio, R = [i"(C 18 0)//(N 2 H+)], 
as a function of mass surface density. Values or R less than 1 
indicates chemically evolved material (see text). The observed R 
drop is consistent with CO being depleted at high densities due 
to the freeze-out. 



ing toward the denser Filament 2 and causing the observed 
blueshift between the C ls O and N2H + lines. The relative ve- 
locity of Filaments 1 and 3 along the line of sight is 3 km s _1 
(Figure|9j. Assuming that the velocity components along the 
plane of the sky have similar magnitudes, this implies a col- 
lision velocity of ~ 5kms _1 . The (slow) shock which fol- 
lowed the impact has caused a density increase in the in- 
teraction region. The (post-shock) density enhancement has 
then contributed to the relatively large average density ob- 
served across G035. 39-00. 33 (Figure fTS)) . Before the collision, 
the filaments appear to be quiescent and coherent in veloc- 
ity, as suggested by the relatively narrow lines and small 
velocity variations observed in the North (Figure[5}, where 
Filament 3 appears spatially separated from Filament 1. It is 
interesting to note that the velocity coherence is maintained 
after collision, as indicated by the presence of velocity gra- 
dients similar to those measured in nearby low-mass dense 
cores (^lkms" 1 pc -1 , see Section [3,2.2l) . 

Note that, in this scenario, shock velocities of about 
5kms _1 in gas of ~10 3 cm -3 appear sufficient to inject 
enough silicon in the gas phase to allow SiO to form and 
be detected across G035. 39-00. 33 (Paper I). However, these 
values of shock velocities and densities are too low for sput- 
tering and grain-grain collisions in magnetized shocks to 
be effective in releasing silicon into the gas phase (e.g. 



Jimenez- S erra et aill2008l . iGusdorf et al.l 



2008, Guillet et al 



201 ll ). It has been discussed in Paper 



I, that in order to reproduce the observed fractional abun- 
dance of SiO via sputtering within C-type shocks, some SiO 
has to be present in the icy mantles (see also Schilke et al. 
1997 and Jimenez-Serra et al. 2008), which need a shock ve- 
locity slightly above 10 kms -1 to be sputtered. Therefore if 
sputtering is indeed the mechanism producing the observed 
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Figure 15. Number density map of G035. 39-00. 33 obtained from 
the observed N2H+ (3 — 2)/(l — 0) line intensity ratio and the use 
of the RADEX code. A constant kinetic temperature of 15 K is 
assumed across the cloud. Some regions at the edge of the filament 
show significantly higher density values than the average value, 
suggesting that the assumption of constant temperature may not 
be valid (see Figure HBT l . Crosses indicate the position of massive 
cores from BT12. 



Figure 16. Kinetic temperature map of G035. 39-00. 33 obtained 
from the observed N2H+ (3 — 2)/(l — 0) line intensity ratio 
and the use of the RADEX code. A constant number density 
of 5xl0 4 cm -3 (the average number density in Figurc[l5} is as- 
sumed across the cloud. The line intensity ratio at the edge of the 
dense filament (Filament 2) is consistent with gas temperatures 
close to 20 K. Crosses indicate the position of massive cores from 
BT12. 



SiO emission, the total velocity in the plane of the sky would 
have to be at least 10 km s _1 . Similarly, in case of grain-grain 
collisions within C-type shocks, the shock velocity has to be 
at least 10 kms -1 to allow dust grains with different elec- 
trical charges to achieve a relative velocity larger than the 
threshold velocity for icy mantle vaporization (6.5 kms , 
Tielens et al. 1994; Caselli et al. 1997). 

Alternatively, the release of SiO molecules mixed in icy 
mantles could be the result of grain mantle vaporization due 
to grain-grain collisions during the process of filament merg- 
ing, without invoking C-type shocks. In this scenario, the 
grain-grain collision velocity can be assumed close to the rel- 
ative velocities of the merging filaments, which is very close 
to the threshold velocity mentioned above, especially taking 
into account the uncertainties in both the observed relative 
velocity and the theoretical threshold velocity value. In fact, 
on the one hand, the inferred relative velocity of ~5 kms -1 



has been obtained assuming that the velocity components 
along the plane of the sky have the same magnitude as 
the component along the line of sight. On the other hand, 
the 6.5 kms -1 value given in Tielens et al. (1994) refers to 
water ice, with a binding energy larger than 5000 K. If icy 
mantles are instead made out of a mixture of H2O and CO 
(as expected in regions with a gas phase dominated by CO 
molecules), binding energy may be significantly lower than 
5000 K and the vaporization threshold velocity will be con- 
sequently lower. A factor of 5 decrease in binding energy 
(appropriate for CO molecules bound on an water ice man- 
tle; Oberg et al. 2005) results in about a factor of 3 decrease 
in the threshold velocity (using the formalism in Tielens et 
al. 1994). Thus, dust grains moving at relative velocities of 
about 5 kms -1 may lose a significant fraction of their icy 
mantles upon collision. A way to test this scenario is to ob- 
serve other species expected to be abundant in the mantles 
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of CO-rich dust grains, in particular CH3OH, which we ex- 
pect to be widespread across the cloud. 



4.2 The Dynamical Evolution of G035.39-00.33 

The high-sensitivity and high-spectral resolution maps in 
both N 2 H+ (1 - 0) and C 18 (1 - 0) have revealed that 
the dense gas of the G035. 39-00. 33 is surrounded by com- 
plex, lower density filamentary structures moving at rela- 
tive velocities of a few kms" 1 . Multiple velocity compo- 
nents have al s o bee n detected in other IRDCs. For instance 
iDevine et al.l (|201 lh found evidence for multiple components 
separated by similar velocities in IRDC G19.30+0.07. In 
the previous sections we have shown that in G035. 39-00. 33, 
these multiple filaments are interacting. 

Given this complexity of the IRDC, the kinematics are 
relatively quiescent. Toward the northern portion of the 
IRDC, where star formation activity is at a minimum, the 
line widths are only a few times the thermal width for a 
15 K gas. Here, C ls O (1 - 0) line widths become close 
to 0.5 kms -1 , similar to those measured in low mass star- 
forming regions (e.g. iFuller fc Myers! 1 19921 ). In the rest of 
the filament, the lines become a factor of 2-4 broader, prob- 
ably due to a combination of filament merging and embed- 
ded (but not yet prominent) star formati on activity (see 
also IDevine et al.l [20lTI . iRagan et al.ll2012l ). Moreover, the 
velocity structure appears coherent, except toward core H6, 
where a lkms -1 pc -1 velocity gradient is measured (prob- 
ably linked to the adjacent zone of interaction). Figure[7] 
shows no large velocity gradients throughout the cloud. 
Hence, the merging of the filaments must be fairly "gen- 
tle", possibly regulated by large-scale magnetic fields. The 
presence of the velocity shift between C ls O and N 2 H + lines 
is telling us that the accumulation of material toward the re- 
gion of interaction of the filaments is still ongoing. Assuming 
a simplified ID geometry, the time to build up the mass of 
Filament 2 from gas flowing in from Filaments 1 and 3 at 
relative velocity w flow is: 



tlio 
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Rf 



Vflo 



0.5 pcy \ 5 kms 1 
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10 



Myr, 



where Rf is the radius of the filament, here normalized to 
0.5 pc, i.e. 36" (see e.g. left panel of Fig[2}, nu,f is the mean 
total hydrogen nuclei number density in Filament 2 and 
fiH,flow is the mean density in the merging regions of Fil- 
aments 1 and 3. The ratio of these mean densities is difficult 
to estimate, especially given the possibility of unresolved 
clumping of the emission of particular gas tracers. From the 
extinction map, Paper III estimated a density contrast of a 
factor of 3 for the case of the "Inner Filament" and its im- 
mediate envelope, which is likely to be a lower limit to the 
ratio used in Equation l4.2l Jimenez-Serra et al. (in prep.) 
find a density contrast of a factor of ~10, based on CO ex- 
citation analysis. The timescale of > 1 Myr that is implied 
for the build-up of Filament 2 by the converging flows of 
Filaments 1 and 3 is longer than the local dynamical time 
in Filament 2 (~ 0.8 Myr), a necessary requirement if it has 
reached virial equilibrium as concluded in Paper III. 



Observer u_ 




Figure 17. Schematic figure illustrating the kinematics of 
G035.39-00.33. 



We have proposed that the merging of filaments has 
led to the local density increase along G035. 39-00. 33 and, 
in particular at the position of core H6. Cluster forma- 
tion due to the collision of filaments has been suggested in 
other star-forming regio ns, for instance in the mass ive star- 
forming region W33A (|Galvan-Madrid et al.ll2010l). in the 
low-m ass star-forming region S erpens ijPuarte-Cabral et alj 
l201ll ). in the L1641-N region (iNakamura et al.ll2012l ). and 
in the Rosette Nebula l|Schneider et alj|2012l ). Here, we have 
been able, using high spectral, high sensitivity data, to show 
that in this particular IRDC, the merging of filaments ap- 
pears to be occurring gently, at a relatively low velocity, and 
that the build-up of material may go on for several million 
years. 



5 SUMMARY & CONCLUSIONS 

We have taken high sensitivity and high spectral resolution 
images of the IRDC G035. 39-00. 33, chosen because of the 
cloud's characteristics suggesting it is in an early stage of 
its evolution. From our analysis we have found the following: 

i) Rather than just the single prominent filament that is 
seen in extinction, we find two additional, morphologically- 
distinct, filaments with velocity components separated by a 
few separated by a few kms -1 . In Filament 2, we find that 
N 2 H + (1 — 0) traces the whole filament and not just the 
cores. There is a local maximum at the position of core H6, 
and this massive core is spatially coincident with the point 
of intersection of the velocity components. We suggest that 
the merging of Filaments 1 and 3 have produced the IRDC 
(Filament 2).. 

ii) Our kinematic study has revealed large-scale velocity 
coherence (comparable line widths at most positions, lack of 
large velocity gradients etc.) across the whole cloud. This is 
indicative of a "gentle" merging of filaments. We suggest that 
the merging of the filaments in the ongoing cloud formation 
process may be resp onsible for the large- s cale S iO (2 — 1) 
emission seen in the \ Jimenez-Serra et al\ \20lA) study. 

iii) We have also witnessed for the first time, a 
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widespread velocity shift (~ 0.2 kms -1 ) between the N2H 4 " 
(1 — 0) and C ls O (1 — 0). Conversely, in low-mass regions, no 
velocity shift is seen between these two tracers. The presence 
of a velocity shift between these tracers over large scales is 
consistent with ongoing merging of filaments, which started 
> 1 Myr ago. Star formation within this cloud may have 
been locally induced by the IRDC formation process (fila- 
ment merging / cloud-could collision). 

iv) Analysis using radiative transfer code, RADEX, 
gives an average Hb number density across the IRDC of 
about 5xl0 4 cm~ 3 and there is indication of a temperature 
gradient perpendicular to the filament. 

v) We find similarities and differences between the stud- 
ied IRDC and nearby (more quiescent) star-forming re- 
gions. Among the similarities we list: (1) material that ex- 
hibits quiescent kinematics, i.e. coherent velocities, espe- 
cially in the northern regions where star formation activ- 
ity is at a minimum (suggesting that the filaments have 
initially physical properties similar to those found in low- 
mass star-forming regions); (2) small velocity gradients 
(~lkms _1 pc _1 ); (3) fractional abundances with respect to 
H2 molecules of N2H + ~4xl0~ 10 . Among the differences 
we found: (1) non-thermally dominated motions traced by 
N2H + , perhaps due to the fact that the dense gas in the 
IRDC has been and continues to be formed from the col- 
lision of supersonically interacting filaments; (2) large-scale 
emission of N2H 4 " (1 — 0) and (3 — 2), implying relatively 
large densities (riH ~ 10 5 cm~ 3 ) extend over large, parsec- 
scale regions in the IRDC. 

vi) This study has also shown the importance of high 
spectral resolution of different gas tracers to unveil the kine- 
matics and correctly interpret line asymmetries. Asymmet- 
ric blue-peaked profiles found across G035. 39-00. 33 are due 
to multiple velocity components along the line of sight. We 
did not find any evidence of large-scale infall motions. In- 
stead, we find evidence of ongoing accumulation of material 
at the intersection of two merging of filaments. Higher angu- 
lar resolution observations are needed to isolate single cen- 
ters of accretion (i.e. the star-forming dense cores analogous 
to the well-studied low-mass cores). 
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APPENDIX A: THE GUIDED GAUSSIAN FIT 

Given the complexity of the data, we have developed a semi- 
automatic procedure in order to fit and interpret the data, 
dubbed the Guided Gaussian Fit (GGF) . The fitting method 
is based on defining multiple velocity components of gas 
along the line of sight. In the case of G035. 39-00. 33, we ob- 
serve three velocity components of differing morphologies 
(see Section 13.2.11 and Figures \S\ & . Based on the emis- 
sion seen in Figure [5] we define three windows in CLASS: 
42-44 kms -1 , 44-46 kms -1 , and 46-48 kms" 1 . Using these 
three windows as a guide, we then fit three Gaussian pro- 
files to the average spectrum of both species. The method 
here for the C ls O (1 — 0) is straightforward in that the three 
components are seen at numerous positions throughout the 
cloud and therefore are evident in the average spectrum. 
This result can be seen in Figure IA1I Here, the average 
spectrum of C ls O (1 — 0) is displayed with three Gaussian 
profiles overlaid. These peak at velocities of, (Filament 1) 
43.411 kms" 1 , (Filament 2) 45.197 kms" 1 , and (Filament 
3) 46.330 kms" 1 . In the N 2 H+ (1 - 0), Filaments 1 and 3 
are only seen at a relatively small number of positions (com- 
pared to Filament 2). Thus, for the Filament 1 and 3 com- 
ponents of N 2 H + we take the average spectrum over those 
positions where these components are detected, rather than 
over the full cloud, to avoid dilution. Centroid velocities, as 
determined by the GGF, for N 2 H+ (1 - 0) of Filaments 1, 2, 
and 3 are 42.985 kms" 1 , 45.582 kms" 1 , and 46.834 kms" 1 , 
respectively. 
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Figure Al. Average spectra of the C 18 (1 — 0) shown in 
blue. Overlaid in green are the three components associated with 
the average spectrum at velocities, 43.411 kms" 1 , 45.197 kms" 1 , 
and 46.330 kms -1 . 



By taking the line width output from the resulting fit, 
we can calculate the velocity dispersion using, 

Av 



&VLSR 



2^/2 ln(2) 



(Al) 



We can now use this as the velocity window within which 
we search for the three velocity components. CLASS is then 
used to calculate the integrated intensity over these three 
windows. The standard detection threshold for Gaussian 
profiles is the 3<r level. In the case of G035. 39-00. 33, we 
see three typically overlapping components along the line 
of sight. We have therefore held the detection limit for any 
individual component at S/N 9. 

Using this as a guide we now know for any given spectra 
(in C ls O (1 — 0)), whether 1, 2 or 3 Gaussians should be fit- 
ted, and which filaments should be fitted. In order to further 
verify this, the less complicated isolated component of N 2 H + 
(1 — 0) is overlaid on the spectra of the C 18 0. Approximately, 
at any given position where N 2 H + (1 — 0) is observed one 
can also expect to see C ls O (1 — 0). Therefore, our "extra- 
guide" to our kinematic study is the N 2 H + (1 — 0) emission. 
If mutliple components are detected in N 2 H + (1 — 0), and 
this is reflected in the spectra of the C ls O, then the same 
components are recorded in both species. If a component 
detected in C 18 (1 - 0) is not detected in N 2 H+ (1 - 0), 
then the C ls O (1 - 0) emission must be checked using the 
GGF intensity and the gas distribution seen in Figure [5] If 
indeed the feature meets these two conditions, it is deduced 
to be a component in its own right and is fitted as such. 
However, if the feature is indistinguishable as a component 
in its own right, either due to blending of two filaments, or 
perhaps depletion features, it is flagged as a blended line, 
and is removed from the analysis. An example of this fitting 
proceedure can be seen in Figure IA2I where by we plot the 
spectra of C ls O (1 - 0) in blue, and N 2 H+ (1 - 0) in red, at 
offset (-7.5, 21), at approximately the position of core H6. 
The results of the guided Gaussian fit are shown in green 
and yellow, respectively. 



APPENDIX B: GAUSSIAN AND HFS FIT 
RESULTS IN SELECTED POSITIONS 
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Figure A2. Spectra at offset (-7.5, 21). C 18 (1 — 0) is shown as the dashed profile, in blue with the results of the GGF shown in green. 
N2H + (1 — 0) is shown as the solid profile, in red with the results of the GGF in yellow. 



Table Bl. Results of the guided Gaussian fitting method at selected offsets. 



Position Filament T mb (K) II (K kms _i ) Vlsr (kms _i ) Av (kms _i ) 



N 2 H+(l-0) 



(-7.5, 21.0) 


Filament 2 
Filament 3 


0.75 
0.40 


(0.04) 
(0.04) 


0.98 
0.49 


(0.06) 
(0.05) 


45.37 
46.73 


(0.03) 
(0.06) 


1.23 
1.16 


(0.06) 
(0.11) 


(5.5, 21.0) 


Filament 2 
Filament 3 


0.93 
0.38 


(0.04) 
(0.04) 


1.11 
0.61 


(0.03) 
(0.02) 


45.32 
46.52 


(0.01) 
(0.02) 


1.12 
1.52 


(0.02) 
(0.06) 


(5.5, -31.0) 


Filament 2 


0.58 


(0.04) 


0.78 


(0.02) 


45.53 


(0.01) 


1.26 


(0.03) 


(5.5, -83.0) 


Filament 1 


0.80 


(0.03) 


1.10 


(0.02) 


45.35 


(0.01) 


1.29 


(0.02) 








C 


i8 0(l- 


-0) 










(-7.5, 21.0) 


Filament 1 
Filament 2 
Filament 3 


1.27 
1.37 
1.24 


(0.11) 
(0.11) 
(0.11) 


1.55 
1.75 
1.67 


(0.07) 
(0.12) 
(0.09) 


43.66 
45.26 
46.76 


(0.02) 
(0.02) 
(0.03) 


1.15 
1.20 
1.27 


(0.05) 
(0.08) 
(0.07) 


(5.5, 21.0) 


Filament 1 
Filament 2 
Filament 3 


0.78 
1.99 
1.13 


(0.10) 
(0.10) 
(0.10) 


1.04 
2.37 
1.47 


(0.05) 
(0.10) 
(0.08) 


43.59 
45.29 
46.72 


(0.03) 
(0.02) 
(0.03) 


1.25 
1.12 
1.23 


(0.07) 
(0.04) 
(0.07) 


(5.5, -31.0) 


Filament 2 


1.65 


(0.13) 


3.46 


(0.06) 


45.55 


(0.02) 


1.97 


(0.04) 


(5.5, -83.0) 


Filament 2 


2.01 


(0.14) 


3.45 


(0.06) 


45.13 


(0.01) 


1.61 


(0.04) 
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Table B2. Results of the HFS fitting method at selected offsets. 



Position 


T a „ t *Tau (K) 


Vls« (kms L ) 


Av ( kms 1 ) 


T 


T 

± ex 






N 2 H+(1 - 


0) 






(-7.5, 21.0) 
(5.5, 21.0) 
(5.5, -31.0) 
(5.5, -83.0) 


7.01 (0.13) 
9.06 (0.13) 
5.14 (0.17) 
7.70 (0.15) 


45.74 (0.01) 
45.63 (0.01) 
45.48 (0.01) 
45.38 (0.01) 


1.94 (0.02) 
1.70 (0.01) 

1.40 (0.02) 

1.41 (0.01) 


3.48 (0.12) 
3.42 (0.10) 
3.07 (0.23) 
4.28 (0.16) 


5.0 (0.1) 
5.6 (0.1) 

4.6 (0.2) 

4.7 (0.1) 






N 2 H+(3- 


2) 






(-4.0, 23.0) 
(9.0, 23.0) 
(9.0, -29.0) 
(9.0, -81.0) 


0.99 (0.06) 
1.41 (0.10) 
0.61 (0.09) 
0.85 (0.07) 


45.83 (0.08) 
45.43 (0.05) 
45.69 (0.09) 
45.36 (0.09) 


2.40 (0.19) 
1.65 (0.16) 
0.81 (0.30) 
2.02 (0.20) 


0.10 
0.10 
0.10 
0.10 





